Function and dynamics of the intrinsically disordered carboxyl terminus of ß2 adrenergic receptor.

Advances in structural biology have provided important mechanistic insights into signaling by the transmembrane core of G-protein coupled receptors (GPCRs); however, much less is known about intrinsically disordered regions such as the carboxyl terminus (CT), which is highly flexible and not visible in GPCR structures. The ß2 adrenergic receptor's (ß2AR) 71 amino acid CT is a substrate for GPCR kinases and binds ß-arrestins to regulate signaling. Here we show that the ß2AR CT directly inhibits basal and agonist-stimulated signaling in cell lines lacking ß-arrestins. Combining single-molecule fluorescence resonance energy transfer (FRET), NMR spectroscopy, and molecular dynamics simulations, we reveal that the negatively charged ß2AR-CT serves as an autoinhibitory factor via interacting with the positively charged cytoplasmic surface of the receptor to limit access to G-proteins. The stability of this interaction is influenced by agonists and allosteric modulators, emphasizing that the CT plays important role in allosterically regulating GPCR activation.

Ultrafast structural changes direct the first molecular events of vision.

Vision is initiated by the rhodopsin family of light-sensitive G protein-coupled receptors (GPCRs)1. A photon is absorbed by the 11-cis retinal chromophore of rhodopsin, which isomerizes within 200 femtoseconds to the all-trans conformation2, thereby initiating the cellular signal transduction processes that ultimately lead to vision. However, the intramolecular mechanism by which the photoactivated retinal induces the activation events inside rhodopsin remains experimentally unclear. Here we use ultrafast time-resolved crystallography at room temperature3 to determine how an isomerized twisted all-trans retinal stores the photon energy that is required to initiate the protein conformational changes associated with the formation of the G protein-binding signalling state. The distorted retinal at a 1-ps time delay after photoactivation has pulled away from half of its numerous interactions with its binding pocket, and the excess of the photon energy is released through an anisotropic protein breathing motion in the direction of the extracellular space. Notably, the very early structural motions in the protein side chains of rhodopsin appear in regions that are involved in later stages of the conserved class A GPCR activation mechanism. Our study sheds light on the earliest stages of vision in vertebrates and points to fundamental aspects of the molecular mechanisms of agonist-mediated GPCR activation.

Impact of membrane lipid polyunsaturation on dopamine D2 receptor ligand binding and signaling.

Increasing evidence supports a relationship between lipid metabolism and mental health. In particular, the biostatus of polyunsaturated fatty acids (PUFAs) correlates with some symptoms of psychiatric disorders, as well as the efficacy of pharmacological treatments. Recent findings highlight a direct association between brain PUFA levels and dopamine transmission, a major neuromodulatory system implicated in the etiology of psychiatric symptoms. However, the mechanisms underlying this relationship are still unknown. Here we demonstrate that membrane enrichment in the n-3 PUFA docosahexaenoic acid (DHA), potentiates ligand binding to the dopamine D2 receptor (D2R), suggesting that DHA acts as an allosteric modulator of this receptor. Molecular dynamics simulations confirm that DHA has a high preference for interaction with the D2R and show that membrane unsaturation selectively enhances the conformational dynamics of the receptor around its second intracellular loop. We find that membrane unsaturation spares G protein activity but potentiates the recruitment of ß-arrestin in cells. Furthermore, in vivo n-3 PUFA deficiency blunts the behavioral effects of two D2R ligands, quinpirole and aripiprazole. These results highlight the importance of membrane unsaturation for D2R activity and provide a putative mechanism for the ability of PUFAs to enhance antipsychotic efficacy.

Chimeric single α-helical domains as rigid fusion protein connections for protein nanotechnology and structural biology.

Chimeric fusion proteins are essential tools for protein nanotechnology. Non-optimized protein-protein connections are usually flexible and therefore unsuitable as structural building blocks. Here we show that the ER/K motif, a single α-helical domain (SAH), can be seamlessly fused to terminal helices of proteins, forming an extended, partially free-standing rigid helix. This enables the connection of two domains at a defined distance and orientation. We designed three constructs termed YFPnano, T4Lnano, and MoStoNano. Analysis of experimentally determined structures and molecular dynamics simulations reveals a certain degree of plasticity in the connections that allows the adaptation to crystal contact opportunities. Our data show that SAHs can be stably integrated into designed structural elements, enabling new possibilities for protein nanotechnology, for example, to improve the exposure of epitopes on nanoparticles (structural vaccinology), to engineer crystal contacts with minimal impact on construct flexibility (for the study of protein dynamics), and to design novel biomaterials.

SCoV2-MD: a database for the dynamics of the SARS-CoV-2 proteome and variant impact predictions.

SCoV2-MD (www.scov2-md.org) is a new online resource that systematically organizes atomistic simulations of the SARS-CoV-2 proteome. The database includes simulations produced by leading groups using molecular dynamics (MD) methods to investigate the structure-dynamics-function relationships of viral proteins. SCoV2-MD cross-references the molecular data with the pandemic evolution by tracking all available variants sequenced during the pandemic and deposited in the GISAID resource. SCoV2-MD enables the interactive analysis of the deposited trajectories through a web interface, which enables users to search by viral protein, isolate, phylogenetic attributes, or specific point mutation. Each mutation can then be analyzed interactively combining static (e.g. a variety of amino acid substitution penalties) and dynamic (time-dependent data derived from the dynamics of the local geometry) scores. Dynamic scores can be computed on the basis of nine non-covalent interaction types, including steric properties, solvent accessibility, hydrogen bonding, and other types of chemical interactions. Where available, experimental data such as antibody escape and change in binding affinities from deep mutational scanning experiments are also made available. All metrics can be combined to build predefined or custom scores to interrogate the impact of evolving variants on protein structure and function.

Unraveling binding mechanism and kinetics of macrocyclic Gαq protein inhibitors.

G proteins represent intracellular switches that transduce signals relayed from G protein-coupled receptors. The structurally related macrocyclic depsipeptides FR900359 (FR) and YM-254890 (YM) are potent, selective inhibitors of the Gαq protein family. We recently discovered that radiolabeled FR and YM display strongly divergent residence times, which translates into significantly longer antiasthmatic effects of FR. The present study is aimed at investigating the molecular basis for this observed disparity. Based on docking studies, we mutated amino acid residues of the Gαq protein predicted to interact with FR or YM, and recombinantly expressed the mutated Gαq proteins in cells in which the native Gαq proteins had been knocked out by CRISPR-Cas9. Both radioligands showed similar association kinetics, and their binding followed a conformational selection mechanism, which was rationalized by molecular dynamics simulation studies. Several mutations of amino acid residues near the putative binding site of the "lipophilic anchors" of FR, especially those predicted to interact with the isopropyl group present in FR but not in YM, led to dramatically accelerated dissociation kinetics. Our data indicate that the long residence time of FR depends on lipophilic interactions within its binding site. The observed structure-kinetic relationships point to a complex binding mechanism of FR, which likely involves snap-lock- or dowel-like conformational changes of either ligand or protein, or both. These experimental data will be useful for the design of compounds with a desired residence time, a parameter that has now been recognized to be of utmost importance in drug development.

Structural basis of the activation of the CC chemokine receptor 5 by a chemokine agonist.

The human CC chemokine receptor 5 (CCR5) is a G protein-coupled receptor (GPCR) that plays a major role in inflammation and is involved in cancer, HIV, and COVID-19. Despite its importance as a drug target, the molecular activation mechanism of CCR5, i.e., how chemokine agonists transduce the activation signal through the receptor, is yet unknown. Here, we report the cryo-EM structure of wild-type CCR5 in an active conformation bound to the chemokine super-agonist [6P4]CCL5 and the heterotrimeric Gi protein. The structure provides the rationale for the sequence-activity relation of agonist and antagonist chemokines. The N terminus of agonist chemokines pushes onto specific structural motifs at the bottom of the orthosteric pocket that activate the canonical GPCR microswitch network. This activation mechanism differs substantially from other CC chemokine receptors that bind chemokines with shorter N termini in a shallow binding mode involving unique sequence signatures and a specialized activation mechanism.

Electron Paramagnetic Resonance Gives Evidence for the Presence of Type 1 Gonadotropin-Releasing Hormone Receptor (GnRH-R) in Subdomains of Lipid Rafts.

This study investigated the effect of type 1 gonadotropin releasing hormone receptor (GnRH-R) localization within lipid rafts on the properties of plasma membrane (PM) nanodomain structure. Confocal microscopy revealed colocalization of PM-localized GnRH-R with GM1-enriched raft-like PM subdomains. Electron paramagnetic resonance spectroscopy (EPR) of a membrane-partitioned spin probe was then used to study PM fluidity of immortalized pituitary gonadotrope cell line αT3-1 and HEK-293 cells stably expressing GnRH-R and compared it with their corresponding controls (αT4 and HEK-293 cells). Computer-assisted interpretation of EPR spectra revealed three modes of spin probe movement reflecting the properties of three types of PM nanodomains. Domains with an intermediate order parameter (domain 2) were the most affected by the presence of the GnRH-Rs, which increased PM ordering (order parameter (S)) and rotational mobility of PM lipids (decreased rotational correlation time (τc)). Depletion of cholesterol by methyl-ß-cyclodextrin (methyl-ß-CD) inhibited agonist-induced GnRH-R internalization and intracellular Ca2+ activity and resulted in an overall reduction in PM order; an observation further supported by molecular dynamics (MD) simulations of model membrane systems. This study provides evidence that GnRH-R PM localization may be related to a subdomain of lipid rafts that has lower PM ordering, suggesting lateral heterogeneity within lipid raft domains.

Tethered agonist exposure in intact adhesion/class B2 GPCRs through intrinsic structural flexibility of the GAIN domain.

Adhesion G protein-coupled receptors (aGPCRs)/family B2 GPCRs execute critical tasks during development and the operation of organs, and their genetic lesions are associated with human disorders, including cancers. Exceptional structural aGPCR features are the presence of a tethered agonist (TA) concealed within a GPCR autoproteolysis-inducing (GAIN) domain and their non-covalent heteromeric two-subunit layout. How the TA is poised for activation while maintaining this delicate receptor architecture is central to conflicting signaling paradigms that either involve or exclude aGPCR heterodimer separation. We investigated this matter in five mammalian aGPCR homologs (ADGRB3, ADGRE2, ADGRE5, ADGRG1, and ADGRL1) and demonstrate that intact aGPCR heterodimers exist at the cell surface, that the core TA region becomes unmasked in the cleaved GAIN domain, and that intra-GAIN domain movements regulate the level of tethered agonist exposure, thereby likely controlling aGPCR activity. Collectively, these findings delineate a unifying mechanism for TA-dependent signaling of intact aGPCRs.

Publisher Correction: GPCRmd uncovers the dynamics of the 3D-GPCRome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

GPCRmd uncovers the dynamics of the 3D-GPCRome.

G-protein-coupled receptors (GPCRs) are involved in numerous physiological processes and are the most frequent targets of approved drugs. The explosion in the number of new three-dimensional (3D) molecular structures of GPCRs (3D-GPCRome) over the last decade has greatly advanced the mechanistic understanding and drug design opportunities for this protein family. Molecular dynamics (MD) simulations have become a widely established technique for exploring the conformational landscape of proteins at an atomic level. However, the analysis and visualization of MD simulations require efficient storage resources and specialized software. Here we present GPCRmd (http://gpcrmd.org/), an online platform that incorporates web-based visualization capabilities as well as a comprehensive and user-friendly analysis toolbox that allows scientists from different disciplines to visualize, analyze and share GPCR MD data. GPCRmd originates from a community-driven effort to create an open, interactive and standardized database of GPCR MD simulations.

HomolWat: a web server tool to incorporate 'homologous' water molecules into GPCR structures.

Internal water molecules play an essential role in the structure and function of membrane proteins including G protein-coupled receptors (GPCRs). However, technical limitations severely influence the number and certainty of observed water molecules in 3D structures. This may compromise the accuracy of further structural studies such as docking calculations or molecular dynamics simulations. Here we present HomolWat, a web application for incorporating water molecules into GPCR structures by using template-based modelling of homologous water molecules obtained from high-resolution structures. While there are various tools available to predict the positions of internal waters using energy-based methods, the approach of borrowing lacking water molecules from homologous GPCR structures makes HomolWat unique. The tool can incorporate water molecules into a protein structure in about a minute with around 85% of water recovery. The web server is freely available at http://lmc.uab.es/homolwat.

DIMERBOW: exploring possible GPCR dimer interfaces.

MOTIVATION: G protein-coupled receptors (GPCRs) can form homo-, heterodimers and larger order oligomers that exert different functions than monomers. The pharmacological potential of such complexes is hampered by the limited information available on the type of complex formed and its quaternary structure. Several GPCR structures in the Protein Data Bank display crystallographic interfaces potentially compatible with physiological interactions. RESULTS: Here, we present DIMERBOW, a database and web application aimed to visually browse the complete repertoire of potential GPCR dimers present in solved structures. The tool is suited to help finding the best possible structural template to model GPCR homomers. AVAILABILITY AND IMPLEMENTATION: DIMERBOW is available at http://lmc.uab.es/dimerbow/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Control of glutamate release by complexes of adenosine and cannabinoid receptors.

BACKGROUND: It has been hypothesized that heteromers of adenosine A2A receptors (A2AR) and cannabinoid CB1 receptors (CB1R) localized in glutamatergic nerve terminals mediate the integration of adenosine and endocannabinoid signaling involved in the modulation of striatal excitatory neurotransmission. Previous studies have demonstrated the existence of A2AR-CB1R heteromers in artificial cell systems. A dependence of A2AR signaling for the Gi protein-mediated CB1R signaling was described as one of its main biochemical characteristics. However, recent studies have questioned the localization of functionally significant A2AR-CB1R heteromers in striatal glutamatergic terminals. RESULTS: Using a peptide-interfering approach combined with biophysical and biochemical techniques in mammalian transfected cells and computational modeling, we could establish a tetrameric quaternary structure of the A2AR-CB1R heterotetramer. This quaternary structure was different to the also tetrameric structure of heteromers of A2AR with adenosine A1 receptors or dopamine D2 receptors, with different heteromeric or homomeric interfaces. The specific quaternary structure of the A2A-CB1R, which depended on intermolecular interactions involving the long C-terminus of the A2AR, determined a significant A2AR and Gs protein-mediated constitutive activation of adenylyl cyclase. Using heteromer-interfering peptides in experiments with striatal glutamatergic terminals, we could then demonstrate the presence of functionally significant A2AR-CB1R heteromers with the same biochemical characteristics of those studied in mammalian transfected cells. First, either an A2AR agonist or an A2AR antagonist allosterically counteracted Gi-mediated CB1R agonist-induced inhibition of depolarization-induced glutamate release. Second, co-application of both an A2AR agonist and an antagonist cancelled each other effects. Finally, a CB1R agonist inhibited glutamate release dependent on a constitutive activation of A2AR by a canonical Gs-Gi antagonistic interaction at the adenylyl cyclase level. CONCLUSIONS: We demonstrate that the well-established cannabinoid-induced inhibition of striatal glutamate release can mostly be explained by a CB1R-mediated counteraction of the A2AR-mediated constitutive activation of adenylyl cyclase in the A2AR-CB1R heteromer.

Conformational Sensors and Domain Swapping Reveal Structural and Functional Differences between ß-Arrestin Isoforms.

Desensitization, signaling, and trafficking of G-protein-coupled receptors (GPCRs) are critically regulated by multifunctional adaptor proteins, ß-arrestins (ßarrs). The two isoforms of ßarrs (ßarr1 and 2) share a high degree of sequence and structural similarity; still, however, they often mediate distinct functional outcomes in the context of GPCR signaling and regulation. A mechanistic basis for such a functional divergence of ßarr isoforms is still lacking. By using a set of complementary approaches, including antibody-fragment-based conformational sensors, we discover structural differences between ßarr1 and 2 upon their interaction with activated and phosphorylated receptors. Interestingly, domain-swapped chimeras of ßarrs display robust complementation in functional assays, thereby linking the structural differences between receptor-bound ßarr1 and 2 with their divergent functional outcomes. Our findings reveal important insights into the ability of ßarr isoforms to drive distinct functional outcomes and underscore the importance of integrating this aspect in the current framework of biased agonism.

Reduced level of docosahexaenoic acid shifts GPCR neuroreceptors to less ordered membrane regions.

G protein-coupled receptors (GPCRs) control cellular signaling and responses. Many of these GPCRs are modulated by cholesterol and polyunsaturated fatty acids (PUFAs) which have been shown to co-exist with saturated lipids in ordered membrane domains. However, the lipid compositions of such domains extracted from the brain cortex tissue of individuals suffering from GPCR-associated neurological disorders show drastically lowered levels of PUFAs. Here, using free energy techniques and multiscale simulations of numerous membrane proteins, we show that the presence of the PUFA DHA helps helical multi-pass proteins such as GPCRs partition into ordered membrane domains. The mechanism is based on hybrid lipids, whose PUFA chains coat the rough protein surface, while the saturated chains face the raft environment, thus minimizing perturbations therein. Our findings suggest that the reduction of GPCR partitioning to their native ordered environments due to PUFA depletion might affect the function of these receptors in numerous neurodegenerative diseases, where the membrane PUFA levels in the brain are decreased. We hope that this work inspires experimental studies on the connection between membrane PUFA levels and GPCR signaling.

Structural Insights into the Process of GPCR-G Protein Complex Formation.

The crystal structure of the ß2-adrenergic receptor (ß2AR) bound to the G protein adenylyl cyclase stimulatory G protein (Gs) captured the complex in a nucleotide-free state (ß2AR-Gsempty). Unfortunately, the ß2AR-Gsempty complex does not provide a clear explanation for G protein coupling specificity. Evidence from several sources suggests the existence of a transient complex between the ß2AR and GDP-bound Gs protein (ß2AR-GsGDP) that may represent an intermediate on the way to the formation of ß2AR-Gsempty and may contribute to coupling specificity. Here we present a structure of the ß2AR in complex with the carboxyl terminal 14 amino acids from Gαs along with the structure of the GDP-bound Gs heterotrimer. These structures provide evidence for an alternate interaction between the ß2AR and Gs that may represent an intermediate that contributes to Gs coupling specificity.

Ligand Binding Mechanisms in Human Cone Visual Pigments.

Vertebrate vision starts with light absorption by visual pigments in rod and cone photoreceptor cells of the retina. Rhodopsin, in rod cells, responds to dim light, whereas three types of cone opsins (red, green, and blue) function under bright light and mediate color vision. Cone opsins regenerate with retinal much faster than rhodopsin, but the molecular mechanism of regeneration is still unclear. Recent advances in the area pinpoint transient intermediate opsin conformations, and a possible secondary retinal-binding site, as determinant factors for regeneration. In this Review, we compile previous and recent findings to discuss possible mechanisms of ligand entry in cone opsins, involving a secondary binding site, which may have relevant functional and evolutionary implications.

Structural basis for the interaction of the beta-secretase with copper.

The ß-secretase (BACE1) features a unique sulfur rich motif (M462xxxC466xxxM470xxxC474xxxC478) in its transmembrane helix (BACE1-TM) which is characteristic for proteins involved in copper ion storage and transport. While this motif has been shown to promote BACE1-TM trimerization and binding of copper ions in vitro, the structural basis for the interaction of copper ions with the BACE1-TM is still not well understood. Using molecular dynamics (MD) simulations, we show that membrane embedded BACE1-TMs adopt a flexible trimeric structure that binds and conducts copper ions through variable coordination. In coarse-grained (CG) MD simulations, the spontaneous assembly of BACE1-TMs trimers results in a right-handed helix packing arrangement. In subsequent atomistic MD simulations the sulfur rich motif defines characteristic copper ion coordination sites along a constricted partially solvated axial pore. Sliding and tilting of BACE1-TMs along smooth A459xxxA463/464xxA467 surfaces, facilitated by a central P472 induced kink, enables copper ions to alternate between different coordination sites, including the prominent C466 and M470. We shed light into the structural arrangement of BACE1-TM trimers and propose a mechanism for copper ion conduction that might also apply to other proteins involved in metal ion transport.